
Calhoun 

iniutuiiaiul AKliiv« ou tfit Nilvdl Poi($ra{jua(« School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2017-06 

Integration of adaptive optics into High 
Energy Laser modeiing and simuiation 

Puent, Donald 

Monterey, California: Naval Postgraduate School 


http://hdl.handle.net/10945/55520 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http://www.nps.edu/ljbrary 


CsMwun is the Neval Postgraduate School's public access distal repository for 
research mate riels and tnstitutjiooal pubftcatiions created by the NPS community. 
Cathouni is named for Professor of Mathematcs Guy K. CatHiuo, NPS's first 
appoinited — and publi^d — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
MontereVr California USA 93943 




NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY, CALIFORNIA 

THESIS 


INTEGRATION OF ADAPTIVE OPTICS INTO HIGH 
ENERGY LASER MODELING AND SIMULATION 

by 

Donald Puent 
June 2017 

Thesis Advisor: Keith Cohn 

Co-Advisor: Joseph Blau 


Approved for public release. Distribution is unlimited. 




THIS PAGE INTENTIONALLY LEFT BLANK 



REPORT DOCUMENTATION PAGE 

Form Approved 0MB No. 0704-0188 

Public reporting burden for this collection of information is estimated to overage 1 hour per response, including the time for reviewing instruction, searching existing 
data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate 
or any other aspect of this collection of information, including suggestions for reducing this burden to Washington headquarters Services, Directorate for information 
Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction 
Project (0704-0188) Washington DC 20503. 

1. AGENCY USE ONLY (LeoveBlonk) 2. REPORTDATE 3. REPORT TYPE AND DATES COVERED 

June 2017 Master’s Thesis (December 2016 to June 2017) 

4. TITLE AND SUBTITLE 

INTEGRATION OF ADAPTIVE OPTICS INTO HIGH ENERGY LASER MODELING AND 
SIMULATION 

5. FUNDING NUMBERS 

6.AUTHOR(S) 

Donald Puent 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey, CA 93943 

8. PERFORMING ORGANIZATION REPORT 

NUMBER 

9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

N/A 

10. SPONSORING / MONITORING 

AGENCY REPORT NUMBER 

11. SUPPLEMENTARY NOTES 

The views expressed in this document are those of the author and do not reflect the official policy or position of the Department of 
Defense or the U.S. Government. IRB Protocol Number: NA. 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Approved for public release. Distribution is unlimited. 

12b. DISTRIBUTION CODE 

13. ABSTRACT (maximum 200 words) 

We investigate the corrective effect of adaptive optics (AO) on High Energy Laser (HEL) beam propagation through atmospheric 
turbulence. The effect of adaptive optics systems was first integrated into a laser propagation scaling code, ANCHOR, that was 
developed within the Direct Energy (DE) research group at the Naval Postgraduate School (NPS). Using a host of plausible input 
parameters, we then compare the estimated target irradiance increase accomplished by various adaptive optics systems. In various 
cases, adaptive optics is shown to substantially increase the effective range of HELs. 

14. SUBJECT TERMS 

directed energy, lasers, adaptive optics 

15. NUMBER OF 

PAGES 71 

16. PRICE CODE 

17. SECURITY CLASSIFICATION 

OF REPORT 

Unclassified 

18. SECURITY CLASSIFICATION 

OF THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 

OF ABSTRACT 

Unclassified 

20. LIMITATION OF 
ABSTRACT 

uu 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. 239-18 


1 




























THIS PAGE INTENTIONALLY LEFT BLANK 


11 



Approved for public release. Distribution is unlimited. 


INTEGRATION OF ADAPTIVE OPTICS INTO HIGH ENERGY LASER 
MODELING AND SIMULATION 


Donald Puent 
Ensign, United States Navy 
B.S., United States Naval Academy, 2016 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN PHYSICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 2017 


Approved by: Keith Cohn 

Thesis Advisor 


Joseph Blau 
Thesis Co-Advisor 


Kevin B. Smith 

Chair, Department of Physics 



THIS PAGE INTENTIONALLY LEFT BLANK 


IV 



Abstract 


We investigate the corrective effect of adaptive optics on High Energy Laser beam propa¬ 
gation through atmospheric turbulence. The effect of adaptive optics systems was first 
integrated into a laser propagation scaling code, ANCHOR, that was developed within the 
Direct Energy research group at the Naval Postgraduate School. Using a host of plausible 
input parameters, we then compare the estimated time-averaged target irradiance increase 
accomplished by various adaptive optics systems. In various cases, adaptive optics is 
shown to substantially increase the effective range of HELs. 
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CHAPTER ONE 


Introduction 


The U.S. Navy is constantly predicting future threats and contemplating new weapon 
technologies to counter them. One of those technologies is directed-energy (DE) weapons. 
Conventional weapons rely on the kinetic energy of projectiles. High energy lasers 
(HELs), one type of directed-energy weapon, work in a fundamentally new way, using 
electromagnetic radiation to damage or destroy enemy assets. 

High energy lasers offer several advantages over conventional weapons, including the 
delivery of energy at light speed, a low cost per shot, virtually unlimited magazines, 
ability to rapidly engage several targets, and ability to operate discreetly. Despite the 
many advantages, DE weapons also have some disadvantages to conventional weapons, 
including a high susceptibility to performance degradation in challenging atmospheric 
conditions. 

The technology of adaptive optics (AO) has proven useful in mitigating the effects of 
the atmosphere on light beam propagation in the field of astronomy [1], where AO is 
often used to improve image resolution. Likewise, AO shows promise in improving 
HEL performance. To better understand how much adaptive optics can improve HEL 
performance and facilitate sound decision making on whether or not to integrate adaptive 
optics with HELs, we simulate the effect of adaptive optics on laser beam irradiance and 
compare that to cases without adaptive optics. Many simulations have been created that 
accomplish this by using numerical diffraction methods. Uniquely, we incorporate adaptive 
optics into a scaling law code to facilitate large parameter studies. 


1 
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CHAPTER TWO 


Overview of Directed Energy 


Directed-energy weapons include all weapons that project electromagnetic radiation. High 
energy lasers are one of the most promising DE weapons and have already undergone 
some field testing. All conventional lasers use the process of stimulated emission to emit a 
coherent light beam. 


2.1 Laser Technology Overview 

The process of stimulated emission was first introduced in 1917 by Albert Einstein, but it 
would take until 1960 before the first laser was built by Theodore H. Maiman. Stimulated 
emission describes the process of a photon of light perturbing an excited electron, causing 
the electron to drop to a lower energy state and emit a coherent photon of the same 
frequency, polarization, and direction as the initial perturbing photon. For this process 
to amplify light, electrons must be energized into the excited state, creating population 
inversion. The material containing the excited electrons is called the gain medium and 
can be made from a variety of materials. A fundamental component of lasers is the optical 
cavity, which is the arrangement of mirrors that surrounds the gain medium and provides 
the feedback loop for stimulated emission to occur [2]. 

Laser designs have used several methods to achieve population inversion. Chemical lasers 
use chemical reactions while gas lasers typical use an electric current or optical pumping to 
excite gases. In weapons research, a common chemical laser is the deuterium fluoride laser, 
which creates fluorine atoms through combustion and combines deuterium and helium 
atoms with the fluorine atoms to create a stable population inversion. A common gas laser 
is the CO 2 laser, which achieves population inversion through burning a hydrocarbon 
and expanding the hot gas through nozzles. CO 2 lasers are used commercially in welding, 
drilling, and cutting [2]. 

Perhaps the most common designs in use today are solid state lasers (SSLs), which use 
solid state materials as the gain medium. Because the applications for solid state lasers are 
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so broad and numerous, the technology has matured rapidly over the past few decades. 
The ubiquity of conventional solid state lasers make them the most ready asset for near- 
term implementation of DE weapons on the battlefield. Fiber lasers are SSLs that use 
optical fibers as the gain medium, making them more rugged and efficient than previously 
developed SSLs. Slab lasers are a particular form of SSLs that use slabs for the gain medium, 
generally defined as large aspect ratio rectangular cross-section materials. Slab lasers offer 
a potential technological advantage over traditional rod lasers in their power-scalability 
[3]. 

A free electron laser (EEL) uses a relativistic electron beam from a particle accelerator as the 
gain medium. The electrons oscillate between a series of alternating permanent magnets, 
emitting light in the process [4]. The emitted light is captured inside of a resonator cavity 
that is similar to the mirror cavity of a conventional laser. Notably, the EEL does not 
require the use of a fixed gain medium with its associated heating and distortion issues, so 
FELs have the potential to be scaled to high output power. Additionally, FELs can be tuned 
to specific wavelengths and produce high quality light beams [4]. 


2.2 Advantages and Disadvantages of Directed-Energy 
Weapons 

Laser weapons are fundamentally different than conventional weapons. As such, there are 
applications for which laser weapons are superior to conventional weapons and others for 
which conventional weapons are better. The delivery of energy through electromagnetic 
radiation, as opposed to the kinetic energy delivered by massive projectiles in conventional 
weapons, results in precise and fast payload delivery, permitting the engagement of quickly 
maneuvering targets. The variable power output of lasers also allow for a graduated 
response, enabling variations of damage that can result in “soft” or “hard” kills. Laser 
weapons are also desirable for their relatively low cost per shot and effectively infinite 
magazines [2]. 

Despite these advantages, directed-energy weapons do have some drawbacks. They are 
large and expensive to build. The propagation of laser light is highly dependent on 
atmospheric conditions and also requires a line-of-sight to the target. Moreover, a finite 
dwell time is required to disable or destroy targets; harder targets lead to dwell times of 
many seconds in some cases [5]. 
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2.3 Brief History of Directed-Energy Weapons 
Development 

The first conference on military applications of lasers occurred in 1963, three years after 
the first laser was created [6], Since then, the U.S. military has invested in various project 
beginning with Project DELTA (CO 2 gas dynamic laser). Navy Advanced Chemical Laser, 
Army Mobile Tactical Unit (CO 2 electric discharge laser), and the Airborne Laser Lab 
(CO 2 gas dynamic laser) throughout the 1970s. Prom 1980-2000s, the Navy developed the 
MIRACL laser, a deuterium fluoride (DP) chemical laser that was the first megawatt class 
laser in the U.S. [7]. MIRACL successfully completed over 150 tests including shooting 
down missiles and drones [2,6]. Perhaps the most well-known DE project was the hydrogen 
fluoride (HP) chemical Space-Based Laser in the 1980s, which was never completed [8]. 
The chemical oxygen iodine (COIL) Airborne Laser (ABL) was a U.S. Air Porce project 
that began in 1996 and successfully destroyed a ballistic missile in 2010 before the project 
was canceled in 2014 [9]. The Tactical High Energy Laser (THEL) was a partnership 
between Israel and the U.S. Army from 1996-2005 that used a DP sub-megawatt laser and 
successfully shot down 28 Katyusha artillery rockets and 5 artillery shells over the course 
of the project [10]. U.S. Army Project ZEUS began in 1996, which uses a solid state laser 
of power ~10 kW to destroy land mines and Improved Explosive Devices (lEDs) and has 
been used successfully on the battlefield in Iraq and Afghanistan [11]. In 2002, the U.S. 
Air Porce began the Advanced Tactical Laser as a follow on to the ABL that, again, used a 
COIL [2]. The Joint High Power Solid State Laser (JHPSSL) project ran from 2003-2010 
and successfully combined 7 slab lasers to generate a laser of overall power -100 kW and 
spawned follow on projects such as the Maritime Laser Demonstrator and the Robust 
Electric Laser Initiative [6]. Most recently, DE weapon projects include the High Energy 
Liquid Laser (2001-present), the Laser Weapons System (LaWS, 2010-present) shown in 
Figure 2.1 that has seen successful testing on the USS Ponce in the Arabian Gulf, and the 
Solid-State Laser Technical Maturation (SSL-TM) project, which has the goal of getting an 
HEL on a ship for testing and evaluation at sea by FY18 [6]. 
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FIGURE 2.1 Picture of the U.S. Navy Laser Weapons System 
(LaWS) onboard the USS Ponce. Source: [12]. 
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CHAPTERTHREE 


Atmospheric Propagation of High 

Energy Lasers 


Light propagation through a vacuum is a well understood phenomena, but various fields 
of study require additional insight on how that propagation changes in the presence of 
an atmosphere. Several effects such as absorption, scattering, thermal blooming, and 
turbulence combine to diminish the irradiance and/or distort the beam as it propagates. 
Beam alteration is a function of wavelength and output power in addition to atmospheric 
conditions. Absorption, scattering, and turbulence are all linear effects, while thermal 
blooming is nonlinear; these effects combine at high output power to create an upper limit 
on the intensity of a light beam at the target. In the case of high energy lasers (HELs), 
the goal of characterizing these effects is to understand their physical limits in certain 
conditions and to develop a mechanism by which it is possible to mitigate them as much 
as possible in order to maximize the intensity of light on targets. We will look at these 
effects separately. 


3.1 Absorption and Scattering 

Absorption and scattering of photons by particles of gases and aerosols are crucial mecha¬ 
nisms that attenuate light. To understand these phenomena, we begin by examining what 
happens when photons interact with particles. 

The analysis begins with a look at the probability that photons will collide with atmospheric 
particles and thereby be absorbed or scattered. Figure 3.1 shows the total effective area 
that is blocked by gas or aerosol particles in a volume V of the atmosphere. This effective 
area is given by the equation 

Aeff = {N(jA)dz 

where Aeff is the total area blocked by all the particles in V, N is the number density of 
particles, u is the effective cross section area of a single particle, A is the cross section area 
of volume V, and dz is the thickness of volume V. To calculate the probability that photons 
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FIGURE 3.1 A volume of air, V, defined by cross section area, A, 
and width, dz, contains particles that absorb and scatter incoming 
light. Each of these processes contribute to the extinction of light. 
Adapted from [13]. 


traveling through volume V will be scattered or absorbed, we take the ratio Aeflf/A, which 
yields 

probability of collision = {Na)dz. 

It should also be noted that the cross-section a depends on the particle as well as the type 
of interaction being considered. 

The number of photons scattered or absorbed in the beam is proportional to the total 
number of photons multiplied by this probability. Since optical power P is proportional to 
the number of photons, the attenuated power dP becomes 

dP 

— = -iNcT)dz 

where the minus sign indicates attenuation. The solution to this equation is the following: 

P{z) = P(0)e-^^^ (3.1) 

where P(0) is the initial power and z is the coordinate along the beam path. This solution 
is known as Beer’s Law, which simply states that light power attenuates exponentially as it 
propagates through the atmosphere. The attenuation or extinction coefficient is defined as 

€ = N<7. 
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FIGURE 3.2 Plot of fractional transmittance (T = P(z)/P(0)) 
through the atmosphere as a function of wavelength. “Atmo¬ 
spheric windows” exist at high T values where the atmosphere is 
transparent. Even within atmospheric windows, there are band- 
widths of high absorption. Source: [13]. 

This characterization of attenuation is not comprehensive, but is scalable for different 
types of gas and aerosol particles. The attenuation coefficient of each type of particle is 
largely independent of the coefficients of other types of particles present in the propagating 
medium, so the total coefficient is simply the sum of coefficients of each type of particle. 
Additionally, there are separate coefficients for absorption and scattering, expressed as 

^ ~ E CTa Pm + Pa 


where a is the absorption contribution, P is the scattering contribution, and the subscripts 
m and a refer to gas molecules and aerosols, respectively. Each of these terms is a function 
of the wavelength since gas and aerosol particles will be optically transparent to some 
wavelengths and opaque at others. 

Figure 3.2 illustrates the attenuation of light through the atmosphere over a broad range 
of wavelengths. The ranges of relative transparency are known as “propagation windows”. 
The bottom chart in Figure 3.2 also illustrates that even within propagation windows, there 
are narrow wavelength bands of high absorption. This makes characterizing the effects of 
absorption and scattering heavily dependent on the wavelength of the light. Furthermore, 


9 



































FIGURE 3.3 The Gaussian irradiance pattern on the left is unal¬ 
tered with beam power shown in white. On the right, an incoming 
cross wind from the left bends the irradiance pattern into a char¬ 
acteristic crescent shape. Adapted from [14]. 


the attenuation coefficient e is often a function of position. This dependence on position 
arises from its dependence on the density of particles present in the beam path, which 
can change significantly with altitude and other factors. With e as a function of position, 
Equation 3.1 is re-written 


P{z) = P(0) exp 



where z is the position along the beam path. 


3.2 Thermal Blooming 

The effects considered thus far have been linear, meaning the beam power at the target 
scales proportionally with the initial power of that beam at the source. In contrast, thermal 
blooming is a nonlinear phenomenon since its effect is dramatically pronounced as beam 
power increases. 

Thermal blooming is a result of the changing physical characteristics of the air as light 
propagates through it. As a light beam’s energy is absorbed by air, the temperature of the 
air increases. This increase in temperature causes the air to expand as hot air is less dense 
than cold air. The air in the beam’s path therefore becomes less dense where the beam 
power is greater. This density differential causes a variation of the index of refraction 
within the beam cross section, with a lower index where the air is warmer, producing 
a diverging lens effect that widens the beam. Cross winds replace heated air along the 
beam path with cooler air. These cross winds cause the beam to bend upwind where the 
index of refraction of the cooler, incoming air is greater. The effect of cross winds on beam 
irradiance is shown in Figure 3.3. 
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To investigate the magnitude of the effect of thermal blooming, consider a beam propagating 
through a uniform atmosphere. The amount of blooming is approximately characterized 
by the dimensionless thermal distortion factor, Nd- Assuming equal weighting along the 
path, Nd is given by 


Nd 


A^flkPnjaT ^ 
PoCp-f^fwind 


(3.2) 


where k = InjX is the wave number, P is the laser power, nr = dn/dT is the rate at which 
index of refraction changes as a function of temperature, a is the absorption coefficient, T 
is the transmittance, po is the density of air, Cp is the heat capacity of air, D is the beam 
diameter, V^ind is the effective wind speed perpendicular to the beam, and R is the total 
slant path length of the beam. The effect of thermal blooming becomes significant when 
Nd > 25 [15]. 


The higher the thermal distortion factor, the more distorted the beam becomes. Note 
the direct proportionality of Nd with P. This characteristic means that the naive action 
of increasing power to compensate for the effects of thermal blooming may actually 
exacerbate the problem. At high Nd, increasing P will more rapidly decrease the irradiance 
on target. The power that produces maximum irradiance on target is known as the critical 
power, Pcrit> and is estimated by the equation 


crit 



/ -1 y/" 

\a(l — m)} 


(3.3) 


where a and m are parameters that describe the intensity profile of the beam. For a Gaussian 
beam, a = 0.0625 and m = 2 [15]. For a typical laser beam and atmospheric conditions. 
Peril is on the order of hundreds of kilowatts. Figure 3.4 shows how increasing the output 
power of a laser will increase the irradiance on target until P = Perit; beyond that level, 
increasing output power will begin to reduce irradiance on target. Note that while the 
calculation for Perit in Equation 3.3 appears to be proportional to the laser power itself, the 
thermal distortion factor is also directly proportional to laser power P, so the laser power 
factor in Pent cancels. 


3.3 Turbulence 

Atmospheric turbulence begins with random temperature fluctuations that originate from 
convection and wind shear and cause air density inhomogenieties. These inhomogeneities 
cause index of refraction variations that bend and distort light. The air pockets of varying 
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Relative Power 


FIGURE 3.4 Plot of irradiance on target as a function of laser 
power. Continuing to increase laser power will decrease irradiance 
on target beyond the value of Pent- The value of Pent is specific 
to the atmospheric conditions in which the laser is operating. 
Adapted from [16]. 


density, called turbulence cells, affect laser propagation by effects known as beam wander 
and scintillation. 


3.3.1 Characterization of Turbulence 

Characterization of turbulence is complicated by nonuniform weather conditions in the at¬ 
mosphere that results in inhomogeneous and anisotropic turbulence cells [17]. Despite this, 
statistical models can approximate small-scale turbulence to be homogenous, isotropic, and 
ergodic, meaning that turbulence on that scale is independent of local position, direction, 
and time, respectively [15]. The upper size-limit to this homogeneous approximation is the 
outer scale of turbulence, Lq. The magnitude of Lq can be about 20% of the height above 
the ground, but is generally not larger than 100 m [17]. These larger wind eddies transfer 
energy to smaller eddies, creating turbulence cells of various size. This cascading effect 
ends in viscous dissipation at the smallest size-scale, whereby kinetic energy transfers into 
thermal energy. The size of the eddies at which viscous dissipation occurs, characterized 
by the variable /q, changes with temperature, pressure, humidity, and wind speed, but 
the eddies are generally on the order of 1 mm near the surface of the Earth to 1 cm in the 
tropopause [17]. Figure 3.5 illustrates this effect. 
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FIGURE 3.5 The cascading effect of energy as it transfers from 
kinetic energy in large wind eddies of size Lq into thermal energy 
through viscous dissipation of smaller eddies with size /q. 


The most ubiquitous parameter describing turbulence is C^, the refractive index structure 
coefficient, and is governed by pressure and temperature differences that create index of 
refraction differences between two points. An intuitive understanding of comes from its 
relation to the structure function, D„(r), which characterizes how much index of refraction 
changes over a displacement of length r. The structure function is defined as 

Dn{r) = {\n{n) - n(n + r)f) (3.4) 

where ri is the starting position, r is the distance from ri over which the structure function 
is being measured, n is the index of refraction, and <|... |2> denotes the averaging over r. 
A statistical model developed by Kolmogorov defines D„ as a function of [18]: 


0 < r «: /o 
/o r Lo 


(3.5) 


A value on the order of 10 m corresponds to weak turbulence, while a value 
on the order of 10“^^ m“^/^ corresponds to strong turbulence [17]. 
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-17 - 16.5 -16 - 15.5 -15 - 14.5 -14 - 13.5 

FIGURE 3.6 The Hufnagel-Valley model roughly characterizes 
values of the refractive index structure coefficient as a function 
of altitude. Near the surface of the Earth, is larger because the 
influence of convection from the ground on turbulence is more 
significant than it is at higher altitudes. Adapted from [19]. 


One of the standard models that characterizes with respect to height is the Hufnagel- 
Valley Turbulence Model [19], shown in Figure 3.6. This model follows the typical trend 
of larger values near the surface of the Earth where the effects of convection from the 
ground heavily influence atmospheric conditions. 


Another common parameter used to describe turbulence over a path is the Fried coherence 
length, ro- The Fried coherence length is the transverse distance over which a beam of light 
will remain coherent and is a function of a specific path through the atmosphere. Coherent 
light here is generally considered to be light that is in phase within a root-mean-square 
difference of 1 radian [15]. The Fried parameter is given by [20] 


and 


ro 


2.905 

6.88 



[plane wave] 


ro 



[spherical wave] 
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where z is the position along the path, R is the total path length, and the additional term 
in the integral for a spherical wave accounts for the reduced beam area affected by the 
turbulence [15]. If a constant C\ is assumed along the beam path for each of these waves, 
the Fried parameters become 


and 


A6/5 

rn = 0.185—;--—— [plane wave] 

J?3/5(C2)3/5 J 


(3.6) 


'*0 


^■^^J?3/5(C2)3/5 


[spherical wave]. 


(3.7) 


Small Fried parameters indicate strong turbulence. Conversely, larger Fried parameters 
indicate weak turbulence. A laser with a wavelength A « 1 pm represented by a spherical 
wave and fired across a range of R = 5 km would experience ro « 1 cm for strong turbulence 
(Cl = 10“^^ m“^/^) and ro » 2 m for weak turbulence (C^ = 10“^^ m“^/^). In the limit of 
weak turbulence where the diameter of the beam director D is less than the Fried parameter, 
the angular divergence of an ideal beam is approximately diffraction-limited at 


e 


A 

D' 


In the limit of strong turbulence where D > vq, the angular resolution becomes atmospherically- 
limited: 

0 - —. (3.8) 

ro 

Therefore, a laser beam coming from an aperture larger than ro will suffer a significant 
reduction in irradiance due to turbulence. 


An additional parameter closely related to ro is the isoplanatic angle, 6o, which is the angle 
subtended by ro over the path of interest. The isoplanatic angle is given by [21] 

- 3/5 

(3.9) 


00 = 


K 

2.905A:^ J Cliz)z^/^dz 


which for a constant value of becomes 

00 = 0.314-. 

R 

The isoplanatic angle is a useful parameter because it defines a region within which 
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turbulence leads to negligible phase differences. In order to correct for the effects of 
atmospheric turbulence on a laser with adaptive optics, the path of the laser beam must 
fall within the isoplanatic angle of the corrected path [21]. 


Finally, the Greenwood frequency is a measure of how quickly a control system must 
correct for turbulence in the atmosphere. Calculating the Greenwood frequency relies 
on Taylor’s frozen cell hypothesis [22], which assumes that turbulence cells are carried 
unchanged by the wind. The Greenwood frequency is calculated by the expression 


/g = 


0.102k‘ 


K 




3/5 


(3.10) 


where k = 2n jX is the wave number, V^ind is the wind speed, and z is the position of the 
beam along the path. Equation 3.10 becomes 

= 0 . 421 ^ 
ro 

for constant V^ind and along the path. It is important to note that Krind is the speed of 
the air in the plane perpendicular to the axis of laser transmission. In the application of 
laser weapons, the relative motion between a laser and a target will manifest in this term. 
For example, a laser propagating along a path of constant C^, an ro of 1 m, and through a 
Vwind of 20 m/s will encounter a Greenwood frequency of fa ~ 10 Hz. 


3.3.2 Effects of Turbulence on Propagation 

The effect that turbulence has on the laser beam depends on the size of the turbulence 
cells. At the upper limit of cell size where the cells are larger than the width of the 
propagating beam, a phenomenon known as “beam wander” occurs that randomly causes 
the entire beam to move about the aim point. At the lower limit of smaller turbulent 
cells, “scintillation” becomes significant, which describes the laser beam breaking up into 
beamlets by the time is reaches the target. In weak turbulence (C^ ~ 10~^^ beam 

wander is the prominent effect, whereas both effects are present in strong turbulence. 

Temperature gradients in the atmosphere result in refractive index inhomogeneities that 
cause a single beam to break up into multiple smaller beamlets through a time-dependent 
effect known as scintillation. Figure 3.7 is a comparison of two irradiance patterns that 
shows this effect. Diffraction plays a significant role in scintillation. A small turbulence 
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FIGURE 3.7 On the left is the simulated irradiance pattern on 
target of a laser affected by an atmosphere with weak turbulence. 
Conversely, the pattern on the right is greatly distorted by a tur¬ 
bulent atmosphere represented by a value that is four orders 
of magnitude larger. Source: [23]. 


cell in the atmosphere will cause beam spreading according to 0 ^ A// where I is the size 
of the turbulence cell and 9 is the diffraction angle. As the turbulence cells become smaller, 
the angle at which light is diffracted from its intended path becomes larger and the final 
irradiance on target becomes more distorted. The maximum diffraction angle occurs for 
the smallest turbulence cells of size /q, where ~ A/Zq. 


Turbulence can be described by its contribution to the resulting time-averaged spot size 
of a laser on target. This concept will be discussed in further detail in Section 5.1, but 
suffice it to say that smaller spot size contribution terms are better than larger terms for 
directed energy applications. The target spot size in a vacuum is adapted from Breaux’s 
formulation and described by [24] 


0^ = 0.5 




(3.11) 


where m' is characteristic of the beam shape as listed in Table 3.1, is a measure of beam 
quality (M^ = 1 represents an ideal beam), A is the wavelength of the laser beam, z is the 
propagation distance, and D is the laser beam director diameter. The spot size contribution 
of turbulence is a function of and is approximated 



where Dg is the effective aperture size appropriate to the beam shape as listed in Table 3.1. 
The total turbulence spot size contribution incorporates the effects of both scintillation 
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TABLE 3.1 Values for m' and as they relate to various beam 
shapes. Adapted from [24], 



m' 


infinite Gaussian beam 

l/n 

2D^ 

truncated Gaussian beam 

0.9166 

0.925D2 

uniform beam 

0.9202 



and beam wander. Scintillation occurs over scale sizes that are small relative to the cross 
section of a laser. Conversely, beam wander occurs over scale sizes that are large relative 
to the cross section of a laser. When considering large-scale effects to be mitigated by a 
fast steering mirror, reduces to the short-term turbulence contribution term a^, which 
is also a function of and is estimated by 

/ 

0.182{ad/M^f(Delrof for(De/ro) < 3 

= 1 (3.12) 

\^{ad/M^f[(De/rof - 1.18(De/ro)^/^] for(De/ro) > 3 

Table 3.1 contains values for m' and Dj as they relate to various beam shapes. The goal 
of adaptive optics is to compensate for the small-scale effect of scintillation in order to 
reduce a^. 
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CHAPTER FOUR 


Adaptive Optics 


Adaptive optics (AO) systems were developed to mitigate the detrimental effects of at¬ 
mospheric turbulence on light beam propagation. AO systems have a variety of different 
designs, but they all address the issue of turbulence and apply similar logic in correcting for 
its effect. We will explore the conditions under which AO systems are useful and provide 
the details of the system components. 


4.1 The Need for Adaptive Optics 

In the previous chapter, we concluded by quantifying the effects of turbulence in terms of 
spot size contribution. The target spot size in a vacuum, a^, is imposed by diffraction and 
is the lower physical limit of total target spot size. The spot size contribution resulting 
from atmospheric turbulence, a^, is in addition to the diffraction variance. If the condition 
flf » is met, the total spot size is dominated by turbulence and AO can be useful in 
improving beam quality. When a] » al the total spot size is diffraction-dominated, so 
AO would not be as helpful a component of the overall HEL system. Consider typical 
values for an ideal HEL modeled as a truncated Gaussian (m' = 0.9166 and = 0.925D^ 
from Table 3.1, = 1) with a beam director diameter D = 50 cm and wavelength A = 

1 pm propagating R = 5 km through an atmosphere of constant intermediate turbulence, 
= 5 X 10~^^ Using Equation 3.7 with these parameters yields Tq « 5 cm. Going 

on to calculate each component of the spot size terms with Equations 3.11 and 3.12 yields 
approximate values of ~ 4.2 x 10“^ m^ and ~ 1.9 x 10“^ m^. These values indicate 
that turbulence is significant for this reasonable case. 

A similar condition used to determine when AO would be useful is the relation between 
the Fried parameter, ro, and the beam director diameter, D. When ro < D, turbulence is 
significant. This condition is also met with the given parameters, so we can see that this 
HEL operating in this plausible environment experiences significant turbulence that can 
potentially see improvement with the use of adaptive optics. 
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FIGURE 4.1 Two sets of images showing the effect of AO taken 
at the Keck Observatory. For each set, the left image is taken 
without AO and the right image is taken after the AO corrections 
are applied. Top: The center of the Milky Way galaxy. Bottom: An 
image of Uranus. Source: [1]. 


4.2 Adaptive Optics Overview 

Before describing particular AO system components and setups, it is useful to outline the 
logic used in developing these systems. We have established how turbulence perturbs a 
wavefront, which disrupts the focus of the beam on target. In astronomy, the target is an 
imaging sensor and the source is an astronomical object, while in the application of laser 
weapons, the target is an enemy asset and the source is the laser. In both applications, it is 
highly desirable to reconstruct the beam quality of the original source before turbulence 
perturbed the wavefront. If we use a deformable mirror in the optics chain to impose 
the opposite of the phase disturbance induced by turbulence, then we can reverse the 
effects of turbulence. This process requires a means of measuring and characterizing the 
phase of the wavefront perturbed by turbulence. Assuming turbulence-induced wavefront 
perturbations can be measured, an AO system also would require a feedback loop to reshape 
the deformable mirror to cancel as much of these perturbations as possible. Adaptive 
optics systems are often closed loop, meaning the effect of the correction applied by the 
deformable mirror is detected by the wavefront sensor and read into future commands to 
the deformable mirror [21]. 

While AO in laser weapons is still a developing field, the technology has been used for 
several decades on telescopes in astronomy to increase image resolution. Two sets of 
images that show the improvement of AO on beam quality are given in Figure 4.1. 
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4.3 Characterizing Wavefront Error 

Opticians classify several types of monochromatic optical aberrations, including spherical 
aberration, astigmatism, coma, field curvature, and distortion [25]. In the application of 
AO, it is necessary to understand these aberrations mathematically. Zernike polynomials 
form a convenient orthonormal basis set defined on a unit circle, making them useful in 
describing an aberrated wavefront. 


The Zernike polynomials, pictured in Figure 4.2, are grouped into even and odd polynomials 
that are defined as 

= Kip) cos{m0) (4.1) 

for even polynomials and 

ZKip^ = Kip) sm(m0) (4.2) 


for odd polynomials where n and m are the radial and azimuthal mode orders respectively, 
p and 0 are polar coordinates on the unit circle, and i?™ is the radial Zernike polynomial 
defined in Equation 4.3 [26]. For these definitions of Zernike polynomials, positive super¬ 
scripts correspond to even Zernike polynomials and negative superscripts correspond to 
odd Zernike polynomials. It is important to note that even when the superscript for odd 
Zernike polynomials is negative, the m value used to calculate the polynomial is positive. 
The mode orders n and m are non-negative integers constrained by the conditions m < n 
and n- m = even. With these constraints, a sequence of ordered index pairs can be listed 
by increasing the value of n and then stepping through each possible value of m. Often, 
the ordered pairs are sequentially assigned the mode order j. A list of the first 10 Zernike 
polynomials is shown in Table 4.1. 


The index, j, is useful for applying Zernike polynomials to characterize atmospheric 
turbulence. Additionally, modified normalization constants are introduced for reasons of 
simplicity in statistical analysis that redefine Zernike polynomials as [27] 


Zj(r, 0) 


Vn + IR^ir)"^ cos(m0) 
^ Vn + li?™(r)V2 sin(m0) 
^KTlRlir) 


j even, m ^ 0 
j odd, 0 ■ 
m = 0 


The radial Zernike polynomials are defined as [27] 


Kir) 


{n-m)l2 

E 

5=0 


i-mn-s)\ 


^n-2s 


s![(n -t m)/2 - s]![(n - m)/2 - s]! 


(4.3) 
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The mode orthogonality can be succinctly written 


where 


/ 


d^rW(r)ZjZy = 5, 


Ojj' 


W(r) = 


/ 

Ijn 

0 


forr < 1 
forr > 1 


If the phase of an electric field is represented as a three-dimensional surface where peaks 
and troughs represent phase advance and retardation, respectively, the Zernike series used 
to describe that surface can be written as 

^) = Xi 

J 

where p = r/R is the normalized radial coordinate, R is the radius of the circle over which 
the series is defined, and aj are coefficients defined as 

R 

aj={l/R^) J (frW{p)(j){p,e)Zj{p,e). 

0 


Tying the mathematical descriptions to a more qualitative understanding of optical aber¬ 
rations, piston is represented by the first Zernike polynomial (/ = 1) and the x and y 
coordinates of tilt are represented by the second (j = 2) and third (J = 3) Zernike polynomi¬ 
als, respectively. Defocus is represented by the radially symmetric fourth (j = 4) Zernike 
polynomial. The higher order effect of astigmatism is represented by the (j = 5) and (j = 6) 
Zernike polynomials while coma is shown by the (j = 7) and (j = 8) polynomials [28]. 

While AO can presumably correct for some number of Zernike polynomials, the remaining 
wavefront is still represented by the residual error, which is conveniently represented with 
Zernike polynomials by 

Aj = ('^7 - E (I%i7 (4-5) 

i=i 

where is the phase variance and J is the number of corrected Zernike polynomials. 
Noll provides the values of residual error up to A 21 . For a wavefront aberrated by turbulence, 
the residual variance is 

Aj X 0.29447“^/^(D/ro)^/^ 
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FIGURE 4.2 Graphical representations of 21 Zernike polynomials 
listed with the modal numbers n and m used to generate them. 
Source: [29]. 


TABLE 4.1 The radial, n, and azimuthal, m, mode order combina¬ 
tions used to produce the first 10 Zernike mode orders, j. Adapted 
from [27]. 


ym 



7-1 


Z“2 

^2 

7-1 


^3 

^3 

j 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


when 7 > 10 [27]. 


4.4 System Components 

Conventional AO systems consist of three basic components: a wavefront sensor, a control 
computer, and a deformable mirror. The basic steps in the corrective process are to measure 
the phase of the beam with the wavefront sensor, compute the necessary deformable mirror 
compensation, and then apply the compensation. 

While the logic underpinning various AO systems is the same, different applications require 
different configurations of the system components. The diagram on the left of Figure 4.3 
shows the most common type of AO in use today, which is the configuration used for 
astronomical imaging [21]. A guide star of some sort emits coherent light from outside of 
the turbulent atmosphere and along a propagation path within the isoplanatic angle of 
the object. The perturbed light is measured by the wavefront sensor and the conjugate of 
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FIGURE 4.3 Left: The most common AO configuration, used for 
astronomical imaging, in use today. Right: The configuration of 
AO in laser weapon systems requires light returned from the target 
to characterize atmospheric turbulence. Sources: [23], [30]. 


the wavefront error is applied by actuators in the deformable mirror. The correction to 
the guide star light beam is then applied to later wavefronts from both the guide star and 
imaging object entering the optics chain. If this loop occurs faster than the characteristic 
turbulence changes {/bw > fc where fsw is the frequency bandwidth of the AO system 
and fa is the Greenwood frequency), image resolution will improve. 

The configuration used in the application of laser weapons has slight variations from that 
shown on the left of Figure 4.3. Laser weapon AO systems differ primarily in their source 
of reference light and the location of the applied correction in the optics chain. While 
imaging systems use a guide star or other source as a reference to measure turbulence, 
laser weapons measure the wavefront of light returning from a target object. It is not 
necessary to use the return from the high average power laser light itself in the wavefront 
compensation measurements. Instead, a collinear pulsed laser with a different wavelength 
and lower average power but higher peak power can be used. The wavefront of the 
reference pulsed laser will be distorted differently than the wavefront of the HEL, but 
the difference can be compensated for by a gain factor in the correction algorithm [21]. 
Moreover, the deformable mirror in an imaging system applies a correction to incoming 
light shortly before it is imaged. Conversely, the deformable mirror in the AO system of a 
laser weapon preforms the phase front of the outgoing laser beam, effectively using the 
turbulent atmosphere to apply a correcting effect to the beam as it reaches its target. The 
AO configuration generally used in laser weapons systems is shown on the right of Figure 
4.3. In addition to partially correcting the abberation induced by atmospheric turbulence, 
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AO systems can passively mitigate static abberation caused by misalignment, imperfection, 
or damaged surfaces in the optical system. 


4.4.1 Wavefront Sensor 

The phase of a light beam is not directly measurable because photo sensors respond to 
irradiance, which is proportional to the square of the complex field. Several methods have 
been developed to measure the phase of a light beam indirectly by measuring tilt and 
curvature [31], [32]. 

The propagation direction of a light beam can be described by a vector normal to the 
surface of its wavefront. Take, for simplicity, an unaberrated wavefront. If this light beam 
were to pass straight through a lens, it would focus to a spot on the optical axis of that lens. 
If, however, this light beam were to pass through the lens at some small angle from the 
optical axis, the spot on the focal plane would move a distance away from the optical axis 
proportional to the light beam’s tilt. The geometry of light beam tilt is shown in Figure 4.4. 



FIGURE 4.4 A tilted wavefront will focus to a spot on the focal 
plane that is a distance 5y away from the optical axis, z. This shift 
is proportional to the tilt, a. Source: [33]. 


By placing a sensor along the focal plane to detect the position of the spot, the tilt of 
the unaberrated wavefront, or average tilt of an aberrated wavefront, can be directly 
measured. One sensor used to measure the spot position is the quadcell, which consists of 
four detectors arranged like quadrants on two-dimensional axes. The difference between 
the average signal in the two left detectors and the average signal in the two right detector 
gives the horizontal position of the spot. Similarly, the difference between the average 
signal in the two upper detectors and the average signal in the two lower detectors gives the 
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FIGURE 4.5 The Hartmann sensor, sometimes called the 
Hartmann-Shack sensor, breaks a beam into constituent beam- 
lets with an array of lenses and measures the tilt of each beamlet 
with a CCD or similar sensor. Source: [34], 


vertical position of the spot. The Hartmann wavefront sensor, shown in Figure 4.5, extends 
this technique of measuring wavefront tilt by breaking a light beam into constituent 
beamlets with an array of lenses. The focused spots are detected by a CCD or an array of 
quadcells and the wavefront tilts of each beamlet collectively form an approximation for 
the wavefront of the overall beam [31]. 

Another method of measuring tilt (or first derivatives of the wavefront) uses shearing 
interferometry. Interferometry makes use of the relation between interference pattern 
intensity and phase. The first experiment characterizing this wave property of light was 
conducted by Young in his famous double-slit experiment [35]. Through the interference 
of light beams, phase information is translated into intensity variations. A common type of 
interferometer used in adaptive optics is a shearing interferometer, which interferes a light 
beam with a displaced replica [36]. The intensity pattern that results from this interference 
is proportional to the displacement, or shearing distance, multiplied by the cosine of 
the slope, or the tilt [21]. Characterizing the full wavefront requires two measurements 
from shearing in two orthogonal directions to measure the slopes in each direction. By 
measuring the interference pattern and shearing distance, the shearing interferometer 
delivers the wavefront tilt. A general setup and the extracted phase information are shown 
in Figure 4.6. 

While Hartmann-Shack sensors and shearing interferometers measure first derivatives, a 
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Interference 



FIGURE 4.6 The shearing interferometer measures the phase dif¬ 
ference between a light beam and a displaced replica of that beam. 
Source: [37], 


wavefront can also be characterized by its curvature. One approach to measuring wavefront 
curvature directly uses an array of lenses, similar to that of a Hartmann-Shack sensor, to 
break a light beam into constituent beamlets. The difference between signals collected from 
two sensors on either side of the focal plane of the lenses correspond to the local curvature 
of the beam [32]. Curvature measurements are especially useful in adaptive optics systems 
that use bimorph deformable mirrors, which produce curvature deformations [21]. 


4.4.2 Deformable Mirror 

In addition to the deformable mirror, most AO systems have a separate tilt mirror in the 
optics chain to remove the overall light beam tilt, which aside from piston is the lowest 
mode of optical aberration. Some adaptive optics systems also contain a secondary mirror 
that can move along the optic axis of the beam to control focus. 

There are generally two kinds of deformable mirrors: segmented and continuous faceplate, 
shown in Figure 4.7. The number of actuators on segmented mirror plates determines 
the allowed motion of each segment. One actuator per mirror adjusts beam translation 
in the piston mode. Three or more actuators per mirror can adjust the tilt mode at each 
segment. Actuators also fall into two categories: displacement actuators and force actuators. 
Actuators commonly push and pull the mirror surface in a direction transverse to the 
mirror surface. Perhaps the most commonly used are displacement actuators that make 
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Continuous Faceplate 



Segmented Faceplate 


FIGURE 4.7 Schematics of two main kinds of deformable mirrors. 
Top: Some coupling exists between the actuators of a continuous 
faceplate mirror, resulting in a more eomplieated mirror response 
to the application of one actuator. Bottom: The segmented face¬ 
plate mirror pictured has one actuator per segment, which only 
provides correction of the piston mode. Source: [38]. 


use of the piezoelectric effect with materials such as lead zirconate titanate (PZT) and 
lead magnesium niobate (PMN) [39] [40]. In the case of a continuous faceplate mirror, 
each actuator will have some influence on the surrounding surface. The total deformation 
caused by a single actuator is modeled by “influence functions” [ 41]. Another closely 
related concept is coupling, which is the height of the mirror at actuators surrounding the 
actuator that applied pressure to the mirror, but is often expressed as a percent. Mirror 
coupling ranges from zero percent for segmented mirrors, a few percent for a thin mirror, 
to 10-20 percent for commonly used DM [21]. All other things being equal, thicker mirrors 
will have more coupling. Low coupling permits higher order aberration correction, but 
requires a larger number of actuators to implement. High coupling results in smoother 
surfaces transitions that can benefit lower order mode correction. Actuators often respond 
differently in different directions of motion, an effect called hys teresis. The position- 
dependence of the mirror on the applied actuator direction is a concern in any AO system 
with actuator materials that exhibit hysteresis. While it can be reduced by calibration, the 
effect of hysteresis ultimately lowers the control bandwidth of the AO system. 
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FIGURE 4.8 Deformable mirror schematic with a square grid of 
64 piezoelectric actuators. Source: [42], 


Other types of deformable mirrors include bimorph and micro-electro-mechanical (MEM) 
deformable mirrors [43] [44]. Bimorph mirrors use two surfaces that expand and contract 
differently in the presence of an applied electric field and deform the curvature of the mirror 
surface. Recently developed MEM deformable mirrors contain hundreds of actuators with 
high control bandwidths and low hysteresis, all of which are ideal parameters for accurate 
reconstruction of higher mode order aberrations. 


4.4.3 Control Computer 

Adaptive optics systems can be described by their relative number of wavefront sensor 
measurements and actuators. A system with fewer wavefront sensor measurements than 
actuators will experience difficulty in determining appropriate actuator commands and is 
called “underdetermined”. If the system has an equal number of both, it is “determined” 
and each wavefront measurement can be directly paired with an actuator. Most systems 
are “overdetermined” meaning there are more signals coming from the wavefront sensor 
than there are actuators to implement the required deformation. In any case, the control 
computer converts electrical signals from the wavefront sensor into signals to the actuators 
of the deformable mirror. Another distinguishing factor of actuators is the speed of the 
response. Faster responses move actuators toward the correct final position more quickly, 
but too fast a response will cause the actuators to overshoot the final position like an 
underdamped oscillator. However, if the response is not fast enough, the AO system 
will perhaps do more harm than good by correcting for the wrong turbulence conditions. 
Generally, an AO system that does not satisfy the condition fsw > fc will not be effective. 
One goal of the control system design is to optimize this AO system response. 

In order to quantify the requirements of a control system, we break down a perturbed 
wavefront into its constituent sine waves using Fourier analysis and analyze the effects 
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of our control system on each of those waveforms. The AO system will effectively alter 
the amplitude and shift the phase of many of these constituent sine waves. The ratio of 
the output amplitude to the input amplitude at each frequency is called the gain. Note 
that in the control system discussion, frequency (tc)) refers to the frequency of a perturbed 
wavefront component, not the frequency of the light beam. If a single input sine wave 
is r{t) = Asin((x)t), its output is c{t) = AM{o))s\n{o)t + <p{t)) where M(co) is the gain and 
^(t) is the frequency-dependent phase shift. Gain is often written in decibels where 
AfdB = 20logM(co). A plot of these values is called a Bode plot as shown in Figure 4.9 [21]. 

While the AO system will eliminate some of the unwanted aberration, there will still 
be some remaining. This remaining aberration is largely a function of the AO system 
bandwidth, which is the range of frequencies over which the system applies a correction 
that sufficiently eliminates those frequencies from the perturbed wavefront. For further 
analysis, it is useful to work in Hz by converting the sinusoidal Fourier components with 
formula f = co/(2n). With the equation [35] 


f(/) = 


we can develop a spectrum of turbulence component frequencies that is also a function of 
the Greenwood frequency, fa- The portion of the remaining aberration that results from 
temporal factors is the residual temporal variance, which can now be calculated as [22] 

OO 

^temp “ / (4.6) 

0 


where H(f, fsw) is a function describing the filter gain with values ranging from 0 to 1. For 
a perfect step filter, H(/, f^w) = 1 when 0 < / < fsw, indicating that those frequencies 
are eliminated by the AO system and will make no contribution to the residual temporal 
variance. Conversely, H(/, fsw) = 0 for fsw < f < indicating a full contribution from 
each of these frequencies to CTjgmp- Solving the integral in Equation 4.6 for this perfect filter 
yields 


^temp 


A 

Jbw , 


\ 5/3 


(4.7) 


\JUVV / 

for higher-order wavefront variance. While real AO systems will not function as perfect 
step filters, this approximation will prove to be sufficient. 
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FIGURE 4.9 Bode plot of amplitude and phase responses of the fre¬ 
quencies of the perturbed wavefront components to a deformable 
mirror correction. Adapted from [45]. 
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CHAPTER FIVE 


HEL Modeling and Simulation 


There are two primary types of simulations in use for high energy laser propagation: 
one using numerical diffraction methods to propagate the the laser beam using the wave 
equation and the other using scaling laws to find target irradiance. Numerical diffraction 
codes are higher fidelity than scaling laws, but can take a prohibitively long amount of 
time for large parameter studies. The quicker, lower fidelity, simulations using scaling 
laws are accomplished by codes such as ANCHOR, which was developed within the NFS 
Direct Energy research group. The faster speed of these scaling codes facilitates exploring 
large, multidimensional design spaces. To study the efficacy of adaptive optics, scaling 
laws accounting for their benefits will be incorporated into ANCHOR, which will then be 
used to estimate improvements using adaptive optics for HEL engagements.. 


5.1 ANCHOR 

The ANCHOR code models the effects of diffraction, turbulence, and jitter on laser beam 
propagation by approximating their effects on target spot size a^, which is then used to 
estimate target irradiance. Each effect contributes a term to the total spot size such that 
[16] 

2 2 . 2.2 /r -i\ 

Q, — Uj 0,^ ( 5 * 1 ) 


where a|, and are the spot size contributions by diffraction. Jitter, and turbulence, 
respectively. The diffraction contribution, following the previous adaptation of Breaux’s 
[24] treatment in Chapter 3, is given in Equation 3.11. This formulation agrees with 
Gebhardt’s [16] work for an infinite Gaussian beam. Eor simplicity, we look at the case of 
a uniform beam such that 


a 


2 

d 



(5.2) 


where characterizes the beam quality of the source, z is the target range, k = InjX 
is the wavenumber, and D is the beam director diameter. The jitter contribution is given 
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by [16]: 




where iO^'^ is the variance of the single-axis jitter angle and assumed to be isotropic such 
that^0^) = The turbulence contribution, , is modeled two different ways depending 
on how much the laser system is designed to mitigate the effect of turbulence. The general 
expression for the short term turbulence contribution is previously given in Equation 3.12 
and restated below for the case of a uniform beam (D = Dg): 


>•■182 ( f )' 


for^ <3 
for ^ > 3 

To 


(5.3) 


where is the diffraction contribution to the spot size, is the beam quality, D is the 
beam director diameter, and ro is the Fried parameter. This turbulence spot size contribution 
models the effect of turbulence without adaptive optics partial correction. The expression 
accounting for partial correction by adaptive optics will be defined as a function of the 
adaptive optics-corrected Strehl ratio in a later section. 


5.2 Modeling Adaptive Optics 

As previously discussed, the atmospheric turbulence perturbs the wavefront of the laser 
beam along the path of propagation. We define a term for the dimensionless phase variance, 
(7^, which is zero for an unabberated wave. While adaptive optics can greatly reduce this 
term on target, physical limitations in each component will apply an imperfect correction, 
leaving a residual phase variance. Specifically, the wavefront sensor will be limited by its 
signal-to-noise ratio. The remaining phase variance from this component is approximated 
by the expression [46] [35] 

^WFS - SNR2 

where SNR is the signal-to-noise ratio of the detector. This approximation holds for sensors 
that measure tilt from light rays coming from a target at infinity. 


The deformable mirror will be unable to perfectly match the detected abberation due to 
having a finite number of actuators, which is captured in the residual wavefront variance 
term [46] [35] [47] 
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where /c is a constant that characterizes the response of the turbulence spectrum to the 
spectrum of the deformable mirror influence functions, is the actuator spacing, and ro is 
the Fried parameter. This approximation is valid for a square arrangement of actuators. 
Previous research has shown that k = 0.34 is applicable to many systems with Gaussian 
influence. 

Finally, the system response will introduce some delay between when an aberration is 
detected and when the appropriate correction is applied. During this time delay, the 
turbulence condition that the laser beam will actually experience will have changed by 
some amount characterized by the Greenwood frequency, fc- The residual variance in this 
case is then a function of the Greenwood frequency and described by Equation 4.7, which 
is restated here: 

~ [fEwl 

The complete treatment of this term is found at the end of the previous chapter. 

These variance terms can be added to give a total residual wavefront variance due to 
turbulence after adaptive optics correction such that [35] 

= <^WFS + ^temp- 

The total residual variance is especially useful in applying the effect of adaptive optics to a 
model of laser propagation because it can be used to approximate the Strehl ratio with the 
Marechal approximation [48]. The Marechal approximation is 

= = (5.4) 

ideal 

The Strehl ratio is the ratio of laser beam intensity on target in a non-ideal case to intensity 
on target in an ideal case. Notably, the contributing terms to this Strehl ratio account only 
for the partially corrected turbulence, not jitter. The ideal case, then, will be the resulting 
intensity with no turbulence or jitter. Importantly, the Marechal approximation only holds 
for small residual variance, where ^2 « 1 [49]. 

These residual phase variance terms were previously used to predict the Strehl ratios for the 
AO system on the La Silla 3.6 meter telescope under various atmospheric conditions [46]. 
These predicted Strehl ratios are listed alongside the observed Strehl ratios in Table 5.1 from 
a study presented by Beckers. Each column is a collection of data taken for one wavelength. 
The Fried parameter is proportional to so each wavelength will correspond to a 
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TABLE 5.1 Strehl ratios collected from the La Silla 3.6 meter tele¬ 
scope adaptive optics system under various atmospheric condi¬ 
tions compared to the Strehl ratios predicted by the Marechal 
approximation. Source: [46]. 


Spectral Band 

Units 

J 

H 

K 

L 

M 

Wavelength 

jjm 

1.2 

1.68 

2.23 

3.87 

4.75 

^0 

cm 

33 

49 

69 
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171 

rs 

cm 

100 

100 

100 

100 

100 

^^ind 

m/sec 

5/15/25 

5/15/25 

5/15/25 

5/15/25 

5/15/25 

To 

msec 

21/7/4 

31/10/6 

43/14/9 

84/28/17 

107/36/21 

fc 

Hz 

6/19/32 

4/13/22 

3/9/15 

2/5/8 

1/4/6 

Td 

msec 

17 

17 

17 

17 

17 

/bw 

Hz 

20 

20 

20 

20 

20 

(jI 

fit 

rad^ 

2.16 

1.12 

0.63 

0.21 

0.14 

2 

^temp 

rad^ 

.15/.95/2.22 

1.2/1.61/2.3 

.67/.94/1.28 

.22/.30/.43 

.15/.20/.28 

^tot 

rad^ 

2.3/3.1/4.4 

1.20/1.61/2.3 

.67/.94/1.28 

.22/.30/.43 

.15/.20/.28 

‘^predicted 

% 

10/4/1 

30/20/10 

51/39/28 

80/74/65 

86/82/76 
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% 

- 

10 

28 

70 

- 


36 





different value in the presence of the same atmosphere (constant C^). The same adaptive 
optics system, characterized by = 100 and fsw = 20 Hz, is used in each case. This data 
neglects the effect of signal to noise ratio in the AO system. Predicted Strehl ratios are then 
calculated for three different wind velocity conditions at each wavelength and compared to 
the observed Strehl ratio. While this study compares values in an astronomical application, 
the agreement between the predicted and observed Strehl ratios makes this model of the 
effect of AO systems a promising one for simulating AO in HELs. 


5.3 Integration of Adaptive Optics into ANCHOR 

The peak intensity of a laser beam after propagating a distance z is 



where P is the initial laser power and e is the extinction coefficient. For the ideal case, we 
return to the assumption of no jitter or turbulence and impose the constraint = 1 such 
that aideal = This gives an ideal peak intensity of 

P 

fideal ~ 2 ^ 

Tta , , 
ideal 

As previously noted, the peak intensity in our non-ideal case only accounts for the case of 
partially corrected turbulence by adaptive optics, neglecting the effect of jitter. The Strehl 
ratio is then 

T 

^ _ ^peak _ “ideal 

^ideal «fdeal + 

To solve for , we rearrange and arrive at the expression 


a 


2 

t 


^fdeald 

S 


(5.5) 


This is the expression for the adaptive optics-corrected turbulence contribution to the total 
spot size. We use Equation 5.2 to estimate and Equation 5.4 to estimate the Strehl 
ratio. Leading to Equation 5.5, we assumed no jitter and an ideal beam. The effects of jitter 
and a non-ideal beam will contribute to the total spot size through the and terms in 
Equation 5.1. 


By following this series of steps, we are able to estimate the effect of adaptive optics on 
target irradiance. For a clearer picture of the impact made by adaptive optics, we investigate 
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target irradiance in each of the two cases: the corrected case with turbulence contribu¬ 
tion estimated by Equation 5.5 and the uncorrected case with turbulence contribution 
estimated by Equation 5.3. The ratio of the target irradiances from each case is the factor 
of performance improvement or degradation from AO contribution. 
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CHAPTER SIX 


Simulations and Analysis 


The effect of adaptive optics was simulated in four scenarios. The first scenario as¬ 
sumed an ideal beam (M^ = l) and no jitter; the second accounted for reduced beam 
quality (M^ = 3); the third accounted for uncorrected jitter (^jitter = 5 prad) — that is, 
the residual platform jitter after tilt correction; and the fourth scenario accounted for 
both effects simultaneously. For the first scenario, we simulate cases for three values 
for control frequency bandwidth (/bw = 10 Hz, 20 Hz, 50 Hz), three values for actuator 
spacing (vg = 0.015 m, 0.03 m, 0.06 m), and three values for beam director diameter 
(D = 0.1 m, 0.3 m, 0.5 m). For the remaining scenarios, we look only at the uncorrected 
and corrected cases for one set of nominal AO parameters {/bw = 20 Hz, = 0.3 m). The 
actuator spacing distances and control frequency bandwidths are typical values obtained 
from the currently available hardware [50] [51]. We also assume that the signal-to-noise 
ratio for the wavefront sensor is sufficiently high so as not to contribute to residual error. 
For the astronomical application of adaptive optics, this was shown to be a plausible 
assumption when the reference beam is sufficiently bright [46]. The laser power is 30 kW, 
so the effect of thermal blooming is negligible and therefore not included in the simula¬ 
tions. Each irradiance plot has three contours along lines of constant peak irradiance at 
lOOMW/m^ (lOkW/cm^), 50MW/m2 (5kW/cm2), and 20MW/m2 (2kW/cm2). These 
peak irradiance values are estimated requirements for hard kills on hard, medium, and soft 
targets, respectively. As a metric of adaptive optics performance, we give approximate kill 
range increases between the tilt-compensated, or uncorrected, cases and adaptive optics 
compensated cases of sea-level and aerial targets. That range increase will be the distance 
the contour moves between the uncorrected and AO compensated cases for each set of AO 
configurations. Another useful metric of AO performance is in the “Improvement Factor” 
plots, which show the factor by which irradiance increases from the uncorrected cases to 
the AO compensated cases for each point in space. 

The adaptive optics systems we simulate are responding to turbulence profiles gener¬ 
ated by a program called LEEDR. The vertical profile that characterizes the operating 
environment of the following laser simulations is shown in Figure 6.1. 
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FIGURE 6.1 The typical vertical C\ profile of Monterey, CA, cre¬ 
ated by LEEDR. This profile is used for the laser simulations in 
this chapter. 


6.1 Ideal Beam Quality, No Jitter Case 

The first scenario is the ideal case where we assume an perfect beam quality (M^ = l) and 
no uncorrected platform jitter. For this case, we investigate the effect of control frequency 
bandwidth, actuator spacing, and beam director size on turbulence compensation. We 
study the ideal beam quality, no jitter case first to isolate the effect of AO on the turbulence 
contribution, , to the total spot size in Equation 5.1. After considering this ideal case, we 
introduce cases with larger diffraction, (by increasing M^), and jitter, (by increasing 
^jitter)j contributions to identify what level of diffraction and jitter contributions begin to 
dominate overall spot size and reduce the benefit of using AO. 


6.1.1 Effect of Increasing Control Frequency Bandwidth 

To characterize the effect of control frequency bandwidth, we choose the nominal beam 
director size of D = 0.3 m. The uncorrected case (no AO) is shown in Figure 6.2. For this 
beam director size, control frequencies of fsw = 10 Hz, 20 Hz, and 50 Hz are evaluated 
and the resulting AO corrected irradiance plots are shown on the left in Figure 6.3. In 
the following discussion, the results for AO correction at each bandwidth frequency are 
compared to the uncorrected case. When fsw = 10 Hz, the kill range of a hard target 
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FIGURE 6.2 Uncompensated (no AO) irradiance pattern for an 
ideal 30 kW laser beam. 


at sea-level increases by approximately 800 m and the soft target kill range increases by 
approximately 1 km. Increasing fsw to 20 Hz can yield a range increase of 1.5 km for a 
hard target and 2.5 km for a soft target. For the case where fsw = 50 Hz, the kill range 
for both target types approximately doubles. For aerial targets (which experience less 
turbulence along the path) at an altitude of 1 km, the kill range of a hard target increases 
by roughly 0.5 km at fsw = 10 Hz, 1 km at fsw = 20 Hz, and 1.5 km at fsw = 50 Hz. For 
soft aerial targets, kill range increases by roughly 2 km at fsw = 10 Hz and to ranges 
well beyond 5 km when few = 20 Hz or 50 Hz. The target irradiance increase gained by 
AO compensation is also shown by the “Improvement Factor” plots that are shown next 
to their corresponding irradiance plots. At a range of 5 km and an altitude of 1 km, AO 
provides improvement factors of approximately 7 and 10 for fgw = 20 Hz and fsw = 50 Hz, 
respectively. At sea-level, those improvement factors increase to approximately 16 and 35. 
Interestingly, when fsw = 10 Hz, AO does not follow the same improvement factor trend 
as the other two cases, where the largest improvement factors are at sea-level. For this AO 
system, target irradiance actually decreases from approximately 0.5 km to sea-level at a 
range of 5 km. We see this results because the Greenwood frequency, fo, is greater than 
10 Hz from sea-level to an altitude of 0.5 km, which means the atmosphere is changing 
faster than the AO can compute and apply the necessary correction. Generally, as AO 
control frequency bandwidth increases, irradiance on target increases. 


6.1.2 Effect of Increasing Actuator Spacing 

For these simulations, we assess the impact of increasing actuator spacing on laser beam 
performance while holding beam director diameter constant at D = 0.3 m. For every other 
set of simulations, the actuators are assumed to be in a 10 x 10 square grid, where the 
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FIGURE 6.3 Peak irradiance plots of the AO corrected laser beam 
and corresponding “Improvement Factor” plots where the beam 
director diameter is held at a constant nominal value of 0.3 m and 
the control bandwidth varies from fsw - Hz to 50 Hz. 
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actuator spacing and beam director diameter scale proportionally such that = D/10. 
In holding the beam director diameter constant for this particular set of simulations, we 
are evaluating actuator arrangements of 20 X 20, 10 X 10, and 5x5 square grids with 
corresponding actuator spacings of = 15 mm, 30 mm, 60 mm. The irradiance plot of the 
uncompensated case is shown in Figure 6.2. 

The AO compensated irradiance plots are shown in Figure 6.4. For hard targets at sea level, 
kill range increases approximately by 1.5 km when = 15 mm and when = 30 mm. 
When Vs = 60 mm, the range increase is about 1 km. For soft sea level targets, the kill range 
increases by approximately 2.7 km, 2.4 km, and 1.5 km when = 15 mm, 30 mm, 60 mm, 
respectively. The effect of AO is less pronounced for an aerial target at an altitude of 
1 km, where hard target kill ranges increase by roughly 1 km and soft target kill ranges 
increase to well beyond 5 km for each simulated value of actuator spacing. In any case, 
larger actuator separation decreases the effect of AO, but appears to have a lesser impact 
on system performance than does control frequency bandwidth. 


6.1.3 Effect of Increasing Beam Director Diameter and Actuator 
Spacing 

For these simulations, we increase beam director diameter size while holding the actuator 
arrangement constant (10 x 10), consequently increasing actuator spacing such that = 
D/10. We hold the control frequency constant at /bw = 20 Hz and simulate the cases 
where D = 0.1 m, 0.3 m, and 0.5 m. The generated irradiance plots are shown in Figure 
6.5, with uncompensated cases shown on the left and AO compensated cases shown on 
the right. When D = 0.1 m, sea-level hard target kill range changes imperceptibly and 
soft target kill range increases by approximately 200 m. For convenience, we restate the 
case where D = 0.3 m and AO correction yields an increase of 1.5 km for hard targets and 
roughly 2.5 km for a soft target. The largest beam director diameter simulated, D = 0.5 m, 
shows hard and soft target kill ranges more than doubling from uncompensated ranges 
of approximately 1 km and 2 km respectively. Hard aerial targets previously discussed 
see imperceptible kill range changes with AO when D = 0.1 m, but large improvements 
of 1 km and 2.5 km when D = 0.3 m and 0.5 m, respectively. Soft aerial targets also see 
imperceptible kill range changes with AO when D = 0.1 m, but increase well beyond 5 km 
for beam director diameters of size D = 0.3 m and D = 0.5 m. 

Notably, the irradiance patterns for the uncompensated cases change with increasing 
beam diameter. This phenomenon is explained by the factors that contribute to target spot 
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FIGURE 6.4 Irradiance plots and “Improvement Factor” plot of AO 
correction by increasing actuator spacing with constant beam di¬ 
rector diameter, D = 0.3 m. Top: 20x20 actuators (15 mm spacing), 
Middle: 10x10 actuators (30 mm spacing). Bottom: 5x5 actuators 
(60 mm spacing) 
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AO Corrected, D = 0.5 m, fniv = 20 Hz 



1 2 3 4 5 

Target Cross Range (km) 


FIGURE 6.5 Peak irradiance plots of the uncorrected (left) and 
AO corrected (right) laser beams where the control bandwidth is 
held at a constant nominal value of 20 Hz and the beam director 
diameter increases with D = 0.1 m (top), D = 0.3 m (middle), and 
D = 0.5 m (bottom). 
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size, which is inversely related to irradiance. In Equation 5.1, we see that three spot size 
factors contribute: and a^. By increasing actuator spacing rs, and also increasing 

beam director diameter D, the AO is effectively decreasing a^, which increases irradiance. 
However, increasing D also decreases the diffraction spot size which also increases 
irradiance. As such, both effects appear in the AO compensated irradiance plots (right), 
but only the effect of decreasing appears in the uncompensated irradiance plots (left). 
The fact that we see such little range increase with AO compensation when D = 0.1 m 
indicates that the total irradiance is dominated by diffraction, and not turbulence, meaning 
that mitigating turbulence with AO will yield only a small benefit for the tested conditions. 
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FIGURE 6.6 Peak irradiance plots for uncorrected and AO com¬ 
pensated laser beams where /bw = 20 Hz and D = 0.3 m. The 
laser beam source is described by = 3 and has no platform jitter. 
Additionally, the “Improvement Factor” plot shows a maximum 
improvement factor of approximately 6. 


6.2 Non-ideal Beam Quality, No Jitter Case 

For the remaining three scenarios, we will consider only the nominal cases of control 
loop bandwidth fsw = 20 Hz and beam director size D = 0.3 m. This scenario evaluates 
the effect of a non-ideal beam (M^ = 3) on turbulence and AO correction. By analyzing 
Equation 5.2, we see that increasing the value of will increase the diffraction spot size 
contribution to the total spot size, relative to the turbulence contribution, a^. 

As shown in Figure 6.6, a sea-level hard target kill range will increase less than 500 m from 
about 1 km to about 1.3 km while a soft target kill range will increase from about 1.8 km to 
about 2.8 km with AO correction. Aerial targets at a 1 km altitude in this environment will 
see kill ranges change negligibly for hard targets and increase by approximately 200 m for 
soft targets. 


6.3 Ideal Beam Quality with Jitter Case 

Our third scenario investigates the effect of uncompensated jitter (^jitter = 5 prad) on 
adaptive optics correction. From Figure 6.7, we can see that a hard target at sea-level can be 
destroyed at approximately 1 km, which increases to approximately 1.5 km with adaptive 
optics correction. Similarly, we can see that a soft target kill range will increase by about 
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FIGURE 6.7 Peak irradiance plots for uncorrected and AO com¬ 
pensated laser beams where= 20 Hz and D = 0.3 m. The laser 
beam source is described by - 1 and has 5 //rad of uncorrected 
platform jitter. A maximum improvement factor is shown to be 
approximately 8. 


1 km from approximately 2 km to 3 km. Looking again at the effect of AO on aerial target 
kill ranges, we see that it will increase by roughly 100 m for a hard target and 500 m for a 
soft target in a case with ^jitter = 5 //rad of platform jitter. 


6.4 Non-ideal Beam Quality with Jitter 

This final scenario accounts for the effects of both a non-deal beam (M^ = 3) and un¬ 
compensated jitter (^jitter = 5 prad) on adaptive optics correction. In this case, both 
and aj contribute to the total spot size, which will tend to mask the contribution and, 
consequently, the effect of AO on overall laser beam irradiance. 

We can see that for a hard target at sea-level, AO provides a kill range increase of a 
approximately 100 m, and for soft targets an increased kill range of about 400 m. The kill 
range increases for an aerial target are approximately 200 m for hard and soft targets. The 
peak irradiance plots and “Improvement Factor” plot for this case are shown in Figure 6.8. 
Evidently, the real effects of reduced beam quality and jitter strongly influence overall 
beam performance, which will render AO less effective in these cases. 
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FIGURE 6.8 Peak irradiance plots for uncorrected and AO com¬ 
pensated laser beams where /bw - 20 Hz and D = 0.3 m. The 
laser beam source is described by = 3 and has 5 //rad of un¬ 
corrected platform jitter. The maximum improvement factor is 
approximately 5. 
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CHAPTER SEVEN 


Conclusions 


One of the major considerations for HEL weapon use is the performance dependence on 
atmospheric conditions. In some cases, adaptive optics can significantly mitigate the effects 
of atmospheric turbulence on HEL performance. In this study, we developed a rapid method 
to simulate the effect of AO on target peak irradiance, enabling large parameter studies. 
The simple cases we examined with nominal AO system parameters show significant 
range improvements. Notably, the effect of adaptive optics is strongest for targets near the 
surface of the Earth where turbulence is stronger, represented by larger values from 
Eigure 6.1. 

The different scenarios that we examined illustrate the significant impact of laser beam 
quality and platform jitter on overall performance. In the ideal scenario, the turbulence 
contribution to the total spot size was relatively large, which resulted in AO correction 
significantly improving the kill ranges of various targets. As beam quality was reduced and 
platform jitter was introduced, we saw a lesser effect of AO on overall target irradiance 
because the turbulence contribution to the total spot size became relatively smaller. As 
such, it is important to understand that AO has the ability to partially correct for the 
diminishing effects of turbulence on laser beam performance, but the performance increase 
is limited to the effect of turbulence. 

Future implementation of high energy laser weapons should consider the corrective capa¬ 
bility of AO in the design process. This simulation can also inform optimization studies that 
vary AO system parameters and evaluate the benefit of the laser performance improvement 
with respect to limiting factors like cost and size of adaptive optics. 
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